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Abstract 

Details of the calculation of the non-forward BFKL kernel at next-to-leading order 
(NLO) are offered. Specifically we show the calculation of the two-gluon production 
contribution. This contribution was the last missing part of the kernel. Together with 
the NLO gluon Regge trajectory, the NLO contribution of one-gluon production and 
the contribution of quark-antiquark production which were found before it defines the 
kernel completely for any colour state in the t-channel, in particular the Pomeron kernel 
presented recently. 
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1 Introduction 



Talking about the BFKL kernel one usually has in mind the kernel of the BFKL equation [1] 
for the case of forward scattering, i.e. for the momentum transfer t = and vacuum quantum 
numbers in the t-channel. However, the BFKL approach is not limited to this particular case 
and, what is more, from the beginning it was developed for arbitrary t and for all possible t- 
channel colour states. Initially it was done in the leading logarithmic approximation (LLA), 
which means summation of terms of the type [a s In s] n (a s is the QCD coupling constant and 
s is the squared c.m.s. energy). This approximation can provide only qualitative results, 
since it does not fix scales, neither for energy nor for transverse momenta determining the 
running coupling constant a s . Therefore calculation of radiative corrections to LLA seems 
to be a daily need. Unfortunately, till now it is not completed, although the forward BFKL 
kernel in the NLO was found already five years ago [2]. 

The problem of the development of the BFKL approach in the next-to leading approxi- 
mation (NLA) is naturally divided into two parts, in compliance with the representation of 
scattering amplitudes in this approach by the convolution of the impact factors of interacting 
particles with the Green's function of two Reggeized gluons in the t-channel. The impact fac- 
tors describing the scattering of particles by the Reggeized gluons contain all the dependence 
on the nature of the particles and are energy independent. All the dependence on energy 
is defined by the universal (i.e. process independent) Green's function, which is determined 
by the BFKL kernel. For a consistent description of scattering amplitudes one needs to 
know the impact factors with the same accuracy as the kernel. Especially interesting is the 
highly virtual photon impact factor, because it can be calculated from the "first principles" 
in perturbative QCD. Unfortunately, this calculation turned out to be a very complicated 
problem, which is not yet solved, although a noticeable progress has been reached here [3]. 
Recently an important step was done finding the solution of a related problem: the NLO 
impact factor for the transition of a virtual photon in a light vector meson was calculated 
in the case of t — and longitudinal polarizations [4]. The NLO impact factors are known 
also at parton level (i.e. for quarks and gluons) [5]. 

The calculation of the NLO BFKL kernel for the non-forward scattering was not com- 
pleted until recently. We remind that for any colour group representation 1Z in the t-channel 
the kernel is given by the the sum of "virtual" and "real" parts [6]. The "virtual" part is 
universal (i.e. it does not depend on TV) and is expressed through the NLO gluon Regge 
trajectory [7]. The "real" part is related to the particle production in Reggeon-Reggeon colli- 
sions and consists of one-gluon, two-gluon and quark- ant iquark contributions. The first con- 
tribution is expressed through the effective Reggeon-Reggeon-gluon NLO vertex [8]. Apart 
from a colour coefficient, it also is universal and known [10]. Each of last two contributions 
is written as a sum of two terms with depending on R coefficients, at that only one of these 
terms enters in the kernel for the antisymmetric colour octet representation 1Z = 8 a (gluon 
channel), whereas the kernel for the colour singlet representation 1Z — 1 (Pomeron channel) 
contains both terms. For the case of quark-antiquark production both these terms are known 
[9]. Instead, only the piece related to the gluon channel was known for the case of two-gluon 
production [10]. Note that for scattering of physical (colourless) particles only the Pomeron 
channel exists. Nevertheless the gluon channel plays an important role. It is caused by the 
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possibility to use this channel for a check of self-consistency, and, finally, for a proof of the 
gluon Reggeization (see Ref. [11] and references therein). 

Thus, the two-gluon production contribution was the only missing piece in the the non- 
forward BFKL kernel. Now it is calculated and the Pomeron kernel is known [12]. Here 
we present the details of the calculation of the two-gluon contributions and the non-forward 
BFKL kernel at NLO for all possible colour states in the t-channel. Since the quark contri- 
bution to the non-forward kernel is known [9] for any TZ, we shall consider in the following 
only the gluon contribution, i.e. we shall work in pure gluo dynamics. 

In the next Section we present the gluon piece of the gluon trajectory, the general form 
of the "real" contribution to the kernel and its part related to one-gluon production. In 
Section 3 we derive the contribution to the kernel from the two-gluon production and define 
the "symmetric" part of this contribution. The colour group relations used in this Section 
are given in Appendix A. The "symmetric" part of the two-gluon contribution is considered 
in Section 4. The three pieces contributing to this part are calculated in Appendices B, C 
and D, respectively. Finally, in Section 5 the non-forward kernel is discussed. 



2 The "virtual" and "one-gluon" parts of the kernel 

As usual, we utilize the Sudakov decomposition of momenta, denoting pi and p 2 the light- 
cone vectors close to the initial particle momenta pa and pn respectively, so that 2pip 2 = 
{Pa + Vb) 2 = s. We use the conventional dimensional regularization with the space-time 
dimension D = 4 + 2e and the normalization adopted in Ref. [6] . The BFKL equation for 
the Mellin transform G„, of the Green's function G is written as 



6; q) = <itt 2 t (D - 2) (* - &) + / S r \^ {n) ((fl ' f; ql G ™ ^ & & ■ (2 ' 1} 



Here qi and q[ = q% — q , {% = l-i-2) are the Reggeon (Reggeized gluon) momenta, q ~ q± is the 
momentum transfer; q 2 ~ q\ = —q 2 = t; the vector sign is used for denoting components 
of momenta transverse to the p\p 2 plane. The BFKL kernel K.^ has the form 



tf*> & q) = [a; (-#) + u (-q{ 2 )} tftfS^ (q, - q 2 ) + K™ q 2 ; q) , (2.2) 

i.e. it is given by the sum of the "virtual" part, determined by the gluon Regge trajectory 
u(t) (actually the trajectory is j(t) = 1 + u>(t)), and the "real" part, related to particle 
production in Reggeon- Reggeon collisions. In the limit e — > [7] the trajectory is given by 



(2.3) 



where to^(t) is the one-loop contribution, whose expression is 

g 2 N c t fd D ' 2 r ^iV c r(l- e )r 2 (e) 
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Here and in the following r' = f — q , N c is the number of colors, T(x) is the Euler function, 
((n) is the Riemann zeta function, (C(2) = n 2 /6) and g is the bare coupling constant, 
concerned with the renormalized coupling in the MS scheme through the relation 
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The "real" part /C^ of the kernel is related to real particle production in Reggeon-Reggeon 
collisions. It can be presented in the form [6] 



where the "non-subtracted" kernel /C^ A is 

Kj»»(a, & a = <W| y E / 14 (*■ *) (?', «)* 5^5n ■ (27) 

Here 7^ is the operator for projection of two-gluon colour states on the representation TZ; 
a, a' and b, b' are Reggeon colour indices; nn is the number of independent states in TZ; 
lib {.0.1)0.2) is the effective vertex for production of the state J in the collision of Reggeons 
with momenta q 1 = (5pi + q±±, g 2 = — ap 2 + q2±', dcfij is the corresponding phase space 
element; the sum is over all possible states J. For a state J consisting of particles with 
momenta fcj, with the total momentum k = qi — q 2 , we have 

tyj = f^(s A - k 2 )(2nf5 D (k - £ h) n ■ (2-8) 

The intermediate parameter sa in Eq. (2.6) must be taken tending to infinity. The second 
term in the R.H.S. of Eq. (2.6) appears only at the NLO and serves for subtraction of the 
large k 2 contribution, in order to avoid a double counting of this region. At the leading 
order (LO) only one-gluon production does contribute, so that k 2 = 0, Eq. (2.7) does not 
depended on sa and gives the kernel in the leading (Born) order: 

)B (a, 02) o) = ? lg q2) ^ (gi ' q ' 2 



G 

_ g 2 N c c n ( o?02 2 + oiol 2 ^\ / 0Q x 

-J&y^K {o,-02? - q ) • (2 ' 9) 

Here 

Tr(v K T d ®T i -) lW\V K \aa') d „ , . 

<* = N c n R = Mv^T^' < 2 ' 10 » 

are the group coefficients, T d are the colour group generators in the adjoint representation, 
T^ b = —if dab) fdab are the group structure constants. The projection operators and the coef- 
ficients en for all possible representations TZ are given, for completeness, in the Appendix A. 
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As it was already noted, the most interesting representations are the colour singlet (Pomeron 
channel, 1Z — 1 ) and the antisymmetric colour octet (gluon channel 1Z — 8 ). Respectively, 
we have for them 

1 



Ci = 1 , 



c 8a = 2 



(2.11) 



Note that for the symmetric colour octet 1Z = 8 S we have c 8s = c 8a , so that at LO the 
BFKL kernels for the symmetric and anti-symmetric octet representations coincide. Since 
the trajectory u(t) is the eigenvalue of K,^ ("bootstrap" of the gluon Reggeization), the 
same is true for K,^^ (trajectory degeneration). However, generally speaking, it does not 
mean that amplitudes with TZ = 8 S in the t-channel have the Reggeized form (with the same 
trajectory but positive signature), because such a form requires impact factors proportional 
to the eigenfunction corresponding to the eigenvalue u(t). It turns out that for parton 
scattering amplitudes in LLA it is just the case. 

Having the one-gluon production vertices at NLO [8], one can easily calculate the one- 
gluon contribution to the kernel with the NLO accuracy. Retaining only terms giving non- 
vanishing contributions in the e — > limit after integration of the kernel over d D ~ 2 k in a 

— # 

neighbourhood of the singular point k — q\ — q\ = 0, we have [10] 

(91, <12] 0) = , ^n_i S I 9 
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(2.12) 



For arbitrary D this part of the kernel can be found in the last of Refs. [8] (see there 
Eq. (4.10)). Note that the exchange qi ^ ql implies also q <-> —q. Since the colour structure 
of the one-gluon production vertex is not changed at NLO, the coefficients c-n here are the 
same as in Eq. (2.9). 

The remarkable properties of the kernel, subsequent from general arguments, are 

} (0, q 2 ; q) = Kf \q u 0; q) = q 2 ; q) = q; q) = (2.13) 

and 

\qu q 2 ; q) = K™{-& ~& ?) = ^\-q 2 , -q ± ; -q) . (2.14) 

Properties (2.13) mean that the kernel turns into zero at zero transverse momenta of the 
Reggeons and appear as a consequence of the gauge invariance. Properties (2.14) are a 
consequence of the crossing invariance and the gluon identity. 

The symmetry properties (2.14) of /C^((fi, q 2 ; q) are evident from Eq. (2.12). Instead 
properties (2.13) are not so evident, but can be easily checked. 
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3 The two-gluon production contribution 



The new states which appear in the sum over J in Eq. (2.7) at NLO are the two-gluon ones. 
For these (and only for these) states the integral in Eq. (2.7) is logarithmically divergent at 
large k 2 . In this region the two-gluon production vertex factorizes into the product of the 
one-gluon vertices (see below), so that the dependence on sa in Eq. (2.6) is cancelled [6]. 

Lorentz- and gauge- invariant representation of the two-gluon production vertex has been 
obtained in Ref. [13]. It has the form 



(T^T*) . . A^(fa, fa) + (T^)., A^(k 2 , fa) 



(3.1) 



where and di are gluon polarization vectors and colour indices, respectively. Note that 
the tensor A^ u (fa,k 2 ) depends not only on fa, fa, as it is explicitly indicated, but on p±,q± 
and p 2 ,q 2 as well. The explicit expression of the tensor is given in Ref. [13]. Its important 
property is the Abelian-type gauge invariance 

KA^fa, M = k 2 A^(k 1: k 2 ) = . (3.2) 

Another important property of the tensor A^fa, k 2 ) is its transformation law under simul- 
taneous exchange (pi, q\, q, i) < > (p 2 , —q 2 , —q,j) , which we will call L < > R exchange: 

A^(fa, k 2 )\ L ^ R = A vll (k 2 , fa) . (3.3) 

This property and the representation (3.1) guarantee that the vertex 7, - 1 2 (qi, q 2 ) is invariant 
with respect to the L < — > R exchange, as it must be. Taking into account the representation 
(2.7), this invariance provides the symmetry of the non-subtracted kernel with respect to 
the exchange q\ <-> — q 2 , q <-> — q . 

We use the Sudakov decomposition for the produced gluon momenta fa and k 2 (k = 
fa + k 2 = qi — q 2 ) in the form (i = 1, 2) 

fa — fliPi + OLip 2 + fa± , saiiPi = -k 2 j_ = fa 2 , Pi = XiP, x 1 +x 2 = l. (3.4) 

We find that it is convenient to use the same light-cone gauge tip 2 = tifa = for both 
gluons, so that we put 

e=e i± — p 2 . (3.5) 

{fap 2 ) 



In this gauge the vertex takes the form 



7g lGa (?i, <fe) = V *elJ 2 ± [( TdlTd2 ) tJ bafifa, h, fa) + (T d ^)^ b^( gi ; fa, fa 
with the tensor 



(3.6) 



1 CX. f-l 

4b<*( qi - t fa, fa) = (gT - P^)(gT - j^M^fa, fa) , (3.7) 

{P2K1) {P2K2) 
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The explicit form of this tensor in terms of the Sudakov variables has been found in Ref. [14]: 



qf± - 2qi±h± 

-x 2 — 



t 



i 



1 / 2 A ± (2x 1 x 2 k ± - A±(x 1 - x 2 ))\ 

¥ [ 2qi±A± + 9ll s ) 



where is the metric tensor in the transverse plane: 

9± - 9 ( P m) ■ { } 

Moreover, the following positions hold: 

A 2 

Aj_ = (x 2 ki - x 1 k 2 )± , k 2 = — , £ = -(xifcai. + £2&i_l) = -A 2 - x\X 2 k\ , 

1 1 

h = (qi-h) 2 = — (x 1 (q 1 -k 1 ) 2 L +x 2 kl ± ) , t 2 = {q\-k 2 ) 2 = —(x 2 (qi-k 2 ) 2 L +x 1 kl ± ) . (3.10) 
X\ x 2 

Note that on account of Eq. (3.7) the tensors b a/3 (qi] k±,k 2 ) and b^ a (q[; k 2 ,ki) contain in the 
denominators "extra" powers of X\ and x 2 , in comparison with "naive" expectations based 
on the consideration of Feynman diagrams. 

We use the notations which coincide with those adopted in Ref. [10]. Note, however, that 
in Ref. [10] two different gauges were used for the two gluons: the gauge (3.5) for the second 
gluon (that with momentum k 2 ), and the gauge similar to the gauge (3.5) with p 2 replaced 
by pi for the first one (that with momentum ki). Therefore our tensors are related to those 
used in Ref. [10] by the equalities 

b al3 (qi; h, k 2 ) = ^(k^ih, k 2 ) , k 2 , h) = ^(A^cy (fc 2 , h) , (3.11) 

where the tensors 



M\k 



n^(k)=gf-2^, (3.12) 

K ± 

with the property 

QT'(k)af(k)=g? , (3.13) 
realize the transformations between the gauges with the gauge fixing vectors pi and p 2 . 

Let us consider the behaviour of the vertex (3.6) in the multi-Regge kinematics, i.e. in 
the limits X\ — > 1 , x 2 — > and X\ — > , x 2 — > 1. The first of them corresponds to the case 
when the first gluon is much closer to the particle A in rapidity space than the second gluon. 
Therefore in this limit the two-gluon production vertex must be factorized as 

7g lga (gi, Q2) = T^gi, gi ~ fci) {gi r^r 7g 3 (gi - fcx, q 2 ) , (3.14) 
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where 7^(?i, q-i) is the one-gluon production vertex. Indeed, at x\ = 1 , x 2 = we have 

£ = -&2_L , ^2^ 2 = — ^2± 5 *1 = (?1 ~ ^l)± > ^2*2 = > (3.15) 



that gives 
so that 

where 
b afS (qi;h,k 
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G1G2 

ij 



^1-^2,^)^ = 0, 

(q 1 ,q 2 )=4g 2 e"J 2 3 ± (t*T*«) . A*, fc 2 )U 1=1 , 



: 2j|zi=l - 



(91 - feOi 

Since in the gauge (3.5) we have 



qi± - T2~h± 
n 1± 



, (qi± - h±) 

— Ki I -= 



qi± - «i± 



^2± 



2_L 



(3.16) 
(3.17) 

. (3.18) 



7y 1 (9i,9i-*i) = -2i7e 1 a L 7j 



?i± — T2 - K i± 
/c 1± 



(3.19) 



we see that the factorization property (3.14) is fulfilled. Moreover, with account of the result 
(3.18), from this property it follows that 



(gi, g - fci; g) IC^ B (g - fci, g 2 ; g) 
(gi - fci) 2 (gi' - fci) 2 



(3.20) 



This equality can be obtained also directly from Eq. (3.18) using the expression (2.9) for the 
LO kernel. 

In the second limit, i.e. x\ = , x 2 = 1 we get 

£ = -k 2 1± , Xl k 2 = -k 2 l± , xJi = k 2 1± , t 2 = (q 1 - k 2 ) 2 ± , 6 a/J (gi; h, k 2 ) = , 



^ Q (<zi;fc2,fci)L=i 



(gi - 



9i 2 ± , 

fc 2_L 



/3 r 



5i± - &2± - 



A- 2 



■A- 



i_l 



and 



We have also 



7tf lGa (gi, 92) = lu 2 (qi, qi ~ k 2 ) ^ ^ k ^ 2 7 Gl (gi - fc 2, g 2 ) • 



(3.21) 
(3.22) 



( ^ 27r jz^i j (6 a/J (gn fo, fci)M?i; fo, fci)) 

_ (g, g - fc 2 ; q) K™ B (fr - fc 2 , & ; g) 

(?i - &) 2 (?V - fc 2 ) 2 



X2=l 



(3.23) 
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With the help of Eqs. (2.7) and (3.6) the two-gluon contribution to the "non-subtracted" 
kernel is presented in the form 

4g A (<?i, & Q) = / (a R F a (k u k 2 ) + b R F b (k u k 2 )) , (3.24) 

where the group coefficients a R and b R are defined as 

Tr (v n (T dl T d2 ) 0(T dl T d2 )*) Tr (v n (T dl T d2 ) (g)(T d2 T dl )*) 
a-R = Tjh- , b R = — , (3.25) 

and the functions F a and F b as 

F a (k u k 2 ) = b^( qi ; k u k 2 )b a(S (q[; k u k 2 ) + k 2 , kJbpM; k 2 , h) , (3.26) 

F b (k u k 2 ) = b^( qi ; k u k 2 )bp a (q[; k 2 , h) + b^fa k 2 , fc 1 )6 a/3 (gi; fc 2 ) . (3.27) 
It is easy to see that a# = cfj, since 

Tr (Vn(T dl T d2 ) (g)(T dl T d2 )*) = Tr (V n (T dl (g T dl *)(T d ' 2 (g T d2 *)) = 

= Tr (p^T dl (g T dl *V n T d2 (g T d2 *) = -^Tr (f^T dl (g T dl *) (f^T^ 2 (g T d2 *) . (3.28) 

This relation is important for the cancellation of the s A -dependence in the kernel (2.6). Due 
to the result (2.11) it gives, in particular, 

a Q = 1 , a 8a = a 8s = J • (3.29) 



For the coefficients b R , using the relation 

^-pd 1 rpd 2 ^ ^Q(T dl T d2 )* — (T dl T d2 ) (g(T rf2 T a!l )* = if dd2di (T dl T d ' 2 ) ^^T d * = — T rf (^)T 

2 



with account of Eq. (2.10) we obtain 



(3.30) 



b R = a R - \c R = Cr(cr-\) , (3.31) 



2 V 2 

i.e. fci = 1/2, whereas for both symmetric and antisymmetric colour octet representations 
the coefficients bR are zero. This is especially important for the antisymmetric case, since 
the vanishing of b$ a is crucial for the gluon Reggeization. Note that the vanishing of bg g 
means that in pure gluodynamics the kernels for both octet representations coincide at the 
NLO as well as at the LO. 



In terms of the variables k i± and Xi the phase space element d(f>oG can be written as 

hcG = W^ 6{1 ~ Xi ~ X2) n^w^ 5D ~ 2{k± ~ kl± ~ h2±)e{SA " k2) ' (3 ' 32) 



Here the identity of the final gluons is taken into account by the factor 1/2!, so that integra- 
tion must be performed over all the phase space; k± = q±± — g 2 _i_ = q{ ± — q 2± and k 2 must 
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be expressed in terms of Xi and ki± (see Eq. (3.10)). We recall that the parameter s A must 
be taken tending to infinity before the limit e — > 0. From the expression (3.8) one can see 
that the tensors b af3 fall down as l/k 2 at k 2 — > oo at fixed Xj. Therefore the integral over 
hi in Eq. (3.24) is well convergent in the ultraviolet region, so that the restrictions imposed 
by the theta-function can be written as 

k 2 

Xi>^. (3.33) 



Let us discuss the properties (2.13) and (2.14) of the two-gluon contribution to the kernel 
(2.2). As for the subtraction term, its properties (2.13), related to gauge invariance, follow 
directly from the corresponding properties of the "Born" kernel (2.9). Properties (2.14) 
are provided by the appropriate symmetries expressed by Eq. (2.9) and the invariance of 
the logarithm and the integration measure d D ~ 2 r / \r 2 r' 2 ) in Eq. (2.6) under the exchanges 
cfi <-> —q[ , r <-> — f', as well as q\ <-> —q 2 , q — q , r <-> — f . 

Turn now to the non-subtracted contribution (3.24). Since its "gauge-invariance" prop- 
erties (2.13) are provided (see Eqs. (3.26) and (3.27)) by the conversion into zero of the 
tensor b aia2 (qi, hi, k 2 ) at qi± = and at g 2 _i_ = (qi — k± — k 2 )± = 0, they can be easily seen 
from the representation (3.8). The symmetry with respect to the exchanges q} <-> — q{ can be 
also easily seen from the representation (3.24), taking into account Eqs. (3.26) and (3.27). 
As it was already discussed, the symmetry relative to the exchanges q\ <-> — q 2 , q <-> — q 
follows from the general representation (2.7) and the invariance of the vertex 7^ lG2 (gi, q 2 ) 
with regard to the L < — >■ R exchange. However, it is not easy to derive this symmetry from 
the representation (3.24) using Eqs. (3.26), (3.27) and the tensor b a p{qi,ki,k 2 ) presented 
in Eq. (3.8). The matter is that the gauge (3.5) breaks the symmetry between p 2 and p±, so 

that a transformation law of this tensor under the L < >■ R exchange has not a simple form 

like that of Eq. (3.3). Moreover, the choice of the variables x,- L also destroys the symmetry 
between p 1 and p 2 . To restore the symmetry one has to introduce the variables defined 
as 

Vi = — , Vi + V2 = 1 , (3.34) 
a 

concerned with Xi by the relation 

yi _ x 2 kjj_ 
y 2 xi kl ± ' 
Using this relation and taking into account that 

^ 2 = (x 2 h -x x k 2 y = (y 2 k 1 - y x k 2 y 
xix 2 y x y 2 

it is easy to see that the phase space element dipcc (3.32) is invariant under the exchange 
Xi <-> yi. Note that in terms of the Sudakov variables the L < — > R exchange means 

Pi <-> , q 1± <-»• -q 2± , q± <-»• -q ± (3.37) 

or, in terms of x,- L and y^ 

Xi^Vi , qi± <-> -<?2± , q± <-> -g± • (3.38) 



(3.35) 
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As it was already mentioned, under this exchange the transformation law for the tensor b a/3 
is not the same as for A^ v of Eq. (3.3), because of the choice of the gauge (3.5). It is not 
difficult to understand that the transformation low must be 

M<?i; K h)\ L ^ R = Q a J (k 1 )Q f3 s (k 2 )bs 7 (qi; k 2 , h) , (3.39) 

where the tensors f2 a/3 (see the definition (3.12)) take into account the gauge change under 
the L <-> R exchange. One can check directly using Eq. (3.8) that the property (3.39) is 
indeed fulfilled. Together with that of Eq. (3.13), this property demonstrates once more 
that the functions F a and F b , defined by Eqs. (3.26) and (3.27) respectively, are invariant 
with respect to the L <-> R exchange. This means that, due to the invariance of the phase 
space element discussed above, the contribution (3.24) is symmetric under the replacement 

qi <-> -92, q <-> -q- 

From Eqs. (2.6), (3.20), (3.23) and (3.26) it follows that the subtraction term can be 
written as 

„4iu2 2 r fj2+2eu ( 2 \ 

I + *■•(*■• >« (iffi) 

= 8g A N*a R J ^^d-i ( x i F a(h, k 2 )\ xi =i + x 2 F a (k 1 , fc 2 )L=i) 

4^2 . rf2+2eu /*.2\ 

+ J2^ I ^^( F ^^)l^- F ^'^)l-=0 ln (^J • ( 3 - 40 ) 

Here the equality a# = c|., the expression (3.32) for the phase space element and the re- 
strictions given by the inequality (3.33) on Xi were taken into account. Note that the second 
integral in the R.H.S of Eq. (3.40) is completely antisymmetric with respect to the substitu- 
tion q\ <-> — q 2 , q <-> —q. Therefore the subtraction term can be obtained by symmetrization 
of the first integral. Consequently, using the definition 

' IU) ] ^hm-fm+Tr^—dm-m)}, (3.41) 



\x(l — x) J + x (1 — x) 

we can write the two-gluon contribution to the kernel (2.6) in the limit Sa — > oo in the form 

KaGbuto,® - J^yJ-i S J dx ] J^yJ-i { xjl^x) ) + ' (3 ' 42) 

where x = x\ and the operator S symmetrizes with respect to the substitution q\ <-> 
~Q2, q ^ —q- It was used here that F b (ki, k 2 )\ xi=0 = F b (ki, k 2 ) | X1=1 = 0, according to 
the definition (3.27) and the properties b a ^(qi, ki, k 2 )\ xi=0 = b^ a (qi, k 2 , fci)| xi= i = . 

Since the coefficient b R is equal to zero for an octet representation, the first term in 
Eq. (3.42) is determined by the two-gluon contribution to the octet kernel, which is already 
calculated [10]. Therefore, our task is to calculate the second term. However, we have 
found that to calculate just this term is not the most convenient way because of the rather 
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complicated form of the convolution F b . Instead of this we find more convenient to calculate 
the "symmetric" contribution 

ft; 9) = J^yJ-^l dx J J^yJ-i [W^r) + ' (3 3) 

where the function 

F s (h, k 2 ) = F a (k 1: k 2 ) + F b (ki, k 2 ) (3.44) 

is given by the convolution 

(b a(3 (qi; h, k 2 ) + b /3a (q 1 ; k 2 , fci))(6 a/3 (gi; h, k 2 ) + bp a {fa k 2 , h)) . (3.45) 

Taking into account that a 8 = 1/4, we can present the two-gluon contribution to the kernel 
for any representation R as 

£gg((Zi> ft; 9) = 4 ( a « - & r)^gg(9i, ft; 5) + <f 2 ; g). (3.46) 



4 The "symmetric" contribution to the kernel 

Calculation of ICqq^i, q 2 ; q) seems to be more convenient since the sum 



b aP (qi,k 1 ,k 2 )+b^ a (q 1 ;k 2 ,k 1 ) 



„2 L,a uP 
9l_L /i l_L fc 2± 

b 2 h 2 



1 Q /3 qf± ~ 2qi±h± x 2 k^ ± q^ ± - x^^ - fci)J 5i±*oj_(9i ~ 
_ o#± x 2 ~ ~ ry^ 

2 ^2 X 2 t 2 k?,j_t 2 

looks simpler than b a/3 (qi, k±, k 2 ) and b /3a (qi,k 2 ,ki) taken separately. Note, however, that 
F s (ki, k 2 ) does not turn into zero at the points x 1 — and a; 2 = 0, in contrast to F b (ki, k 2 ), 
so that the prescription (3.41) is necessary in Eq. (3.43) to remove the "Regge divergencies", 
i.e. the divergencies at x = X\ = and x — 1. We shall use the decompositions 

& a/3 (<7i; fci, fc 2 ) + 6^(gi; k 2 , h) = bf + bf + bf , 
& a/3 W; h, k 2 ) + 6^^; fc 2 , h) = C 3 + + b' 2 a ? , (4.2) 

with 

h a/3 _ Qi±k?±k 2± , s 

°o - p p ' l 4 -^ 

"'1±"'2± 



2 t x Xi^i fc? , ti 



qf ± - 2q 1± k 1± x 2 k? ± q^ ± - Xiq^fai - h) p ± qf±kf_ L (q 1 -k 1 Y 



11 



and 



= 1 g fx 1 q * ± - 2qi±k2L - x i^ fc 2± - a?2(gi - fc2)igfi _ gi 2 ±(gi - k 2 )V4i_ ^ ^ 

2 ^2 X 2 t 2 k 2 ±t 2 

The tensor 6^ does not depend on x at all; as far as b°^ 2 is concerned, it is easy to obtain 



a/3 1 



\x=0 



ft l±gl_L ^a/3 



k 2 



, w 2 | x =o 



(ft - h)l 



a/3 



X=l /■ . I. \9 7.2 / . 7. \2 ' °: 



*|±(gi - h)l 



(gi - fci): 



kf±(qi - fci) 



a/3 1 

2 |x=l 



/V 2 



(4.6) 



(4.7) 



Note also that at fixed x the tensors have infrared singularities. The singularities of 6q 
are evident from Eq. (4.3). As for fr"' 3 with i — 1, 2, they are singular at where 



a/3 1 



|fci±-»0 



/V 2 



°2 |fe 2 x^0 



gl±^2± 
/V 2 



Accordingly to the composition (4.2), we present F s (ki, k 2 ) in the form 

F 8 (ki, k 2 ) = (A + A 1 + A 2 + A3) + (A* <-> k 2 ) . 



(4.8) 



(4.9) 



Note that, since the total convolution (4.9) includes the terms obtained by the substitution 
{k\ <-> k 2 ) from Aj, the choice of Aj is not unique. We get (in the following in this Section 
only transverse momenta are used and we omit the _L sign; pay attention, however, that the 
Minkowski metric is used) 

-i ,;>. ,/2 

(4.10) 



Ac 



2°0 °0aj3 



g?gf 



Ai = bfb> 0a „ + bfb' la 



Qi 
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k\k 2 t\ 



*/3 



^{hhM - 2q 1 k 1 ) + -kl( qi k 2 ) + (ql - (qikMhiqi ~ h)) 
2 x\ 



T ct/3? / 
b l »la/: 



2*i*i 



(D-2) 



+ (gi <-»• g'i) , 
(4.11) 



x 2 (g 2 - 2gifci)(gf - 2q[ h) + x 2 (q[ 2 - 2q' 1 k 1 ) 



x ( (9i*i)(- + T2) ~ 2gi ) + % 2 (gigi ) + (gigi )((gi " k.M - h)) - 2^(k 1 q 1 )(q[(q l - h)) 



x 2 



"xi k{ 



x\ 



Xi 



/2 

+% (2^k 2 1 (q 1 (q' 1 - h)) + (ql - 2(g 1 fc 1 ))((« Zl - hM - h)) 
Ki \ X\ 



(q t <- q[) , (4.12) 



As for A 3 , its definition is not so simple. Actually A 3 is constructed from b^b^^. Note 
that this convolution is invariant with respect to the simultaneous substitution ki <-> k 2 
and gi «-> q[. As for the two terms in bi b' 2a g, which are obtained from each other by this 
substitution, we have taken one of them in an unchanged form, whereas in the other we have 
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performed the substitution k\ «-> k 2 , so that it can be obtained from the first term by the 
substitution qi <-> q[. Of course, this procedure is not unique. We define A3 as 



2tit f 2 



(D-2) 



x 1 x 2 (q 2 1 - 2q 1 k 1 )(q' 2 - 2q[ k 2 ) + (q[ 2 - 2q[k 2 ) [ (qih)(l + ^) - 2x x q\ 



k\ 



2xo 



HQihMh) + ( qi (q[ - k 2 )M( gi - h)) - -LfafaikM - k 2 )) 

X\ 

, 2 x 2 (q x k 2 ){c( x h - hjky ) - x x (q x k 2 - k x k 2 )(q x 4\ ~ qifa) 



xik 2 

+ T^h^ qik2 ~ ^ 2 )(<?i h - k 2 ki) 

Kl K 2 



+ (Qi <-»■ q[) ■ 



From Eqs. (3.43) and (4.9) it follows that we can present /C^(gi, q 2 ; q) as 

) S(J + J 1 + J 2 + J 3 ) , 



where 



Ji 



dx 



(2tt) d - 1 (4tt) 2 + 

d 2+2 "h (A, + (k 



(4.13) 



(4.14) 



(4.15) 



7r 1+e r(l -e) V x(l-x) J + 

Since A does not depend on x, the integral Jq evidently is equal to zero according to the 
definition (3.41). The integrals Ji with i =^ will be discussed below. Here we note that 
^-gg(<71' 92) Q) * s finite in the limit e — > 0. Indeed, let us consider the terms in Ai having 
non-integrable infrared singularities in the limit e — * . These terms can be easily obtained 
using Eqs. (4.3) and (4.8), as well as the explicit expressions (4.11)-(4.13): 



A l 



q[ 2 (qik 



sing 



2 , 2 • (ft ■ - <h) ■ A 2 



h 2 h 



sing 



Mi) , / m 4 



sing 



. (4.16) 



We see that the infrared singular parts of Ai do not depend on x. On the other hand 
for Aj| x=0 and i4j| x=1 we find (it can be done using Eqs. (4.3), (4.6) and (4.7) as well as 
Eqs. (4.11)-(4.13) ) 



A, 



x=0 



21.2 + fe^ft') ? 



Ai 



x=l 



ft 
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k\k\{q x - h) 



>-2\x=l — 



A 



^4-2 U=0 

(ft - k i)(Q[ ~ fci) 
2(?i - h) 2 (q'i - 

(kM-h)) 



Ml - (Qih))(k 2 (qi - fci)) + (ft <-»• q[) , 

- -i^r + (qi » ft) , 



(4.17) 



J2 



Mi ) + |r (ft 2 ~ 2(*i?i 



+ (ft 



32=0 



(ftVi) - q'i 



, 2 (qi k 2) 



+ (qi <-»• ftO> ^U=i = . 



(4.18) 



(4.19) 



Comparing these expressions with the results (4.16), we see that they do not contain new 
(i.e. different from Ai\ sing ) non-integrable infrared singularities in the limit e — > 0. Therefore 
such singularities are absent in (Ai/[x(l —x)]) + . 
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4.1 Calculation of J\ 



In order to calculate the integral J\ it is suitable to first integrate over x and then over 
k\. Using the invariance of the integration measure with respect to the exchange k\ <-> k 2 , 
after the first integration (here and below we omit in integrands terms giving zero after the 
subsequent integration) we obtain 

^ _ q[ 2 f d 2+2 *h 1 f( qi - h) 2 \ 



Tr^^l-e) k 2 \ k 2 J 

x ( % ^ - t,)» - —) + (9i " ' (4 ' 20) 

Note that the singularities of separate terms in the integrand at ki = 0, k 2 = and 
<?! — fci = cancel each other. 

Details of the calculation of this integral are given in Appendix B. At arbitrary e we find 



n' 2 r 1 r 1 
Ji = ^~ dx dyy*- 1 
2 Jo Jo 



[l-e)qlql x 



+ 



(-(fc 2 (l -x) + q$x)(l -y)- q 2 x{l - x)y) e VI - x 

I _ fc 2 )(l _ x )(l _ y (l _ x )) _ g2 (a . _ y(l _ 3.2)) \ 



-(/c 2 (l - x) + -y) - - x)y) £ x(l -x)J + 



+ (ft <- ft') . (4.21) 



The integral cannot be expressed in terms of elementary functions not only at arbitrary D, 
but even in the limit e — > 0. In this limit we have 



* - f r^f-^V ,^, + (I) >n ( j) j + (. - ,) , 

(4.22) 
where 

T/ , . r 1 , /a(l — rr) + faA 

I(a,b,c)= — r : -. r In [ ■ (4.23) 

jo a(l — x) + ox — cx(l — x) \ cx(l — x) J 

Note that the integral I(a, b, c) is invariant with respect to any permutation of its arguments, 
as it can be seen from the representation [15] 

f 1 f 1 f 1 dx 1 dx 2 dx 3 5(l - x 1 - x 2 - x 3 ) 

I(a,b,c)= / / / - -. (4.24) 

Jo Jo Jo {axi + bx 2 + cx 3 ){xix 2 + x\x 3 + x 2 x 3 ) 

In particular, I(k 2 ,q 2 ,qf) does not change under the substitution ft <-> — g 2 - 



4.2 Calculation of J" 2 



The order of integration used for the calculation of Ji (first over x and then over ki) is 
suitable for the calculation of J 2 as well. Details of the integration are given in Appendix 
C. The result of the integration over x can be presented as 



J2 



d 2+2 *h 



7r 1+e r(l -e) 



(ft - h) 2 - k\ 



9? 




- fci)(gj ~ fci) 
(qi - h) 2 
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/ gf_ _ q 2 \ _ (gig£) 



2 % 



2 



1 + e 



(5i - fci) 5 



i , f(qi-k 



\2' 



^ + (ft - 9i) • ( 4 - 25 ) 



At arbitrary e the integration in Eq. (4.25) gives 



where 



r / / \ _ / ^fel (91 - fcl)(gj - fcl) ]n Agl - ^l) 2 (gl ~ hf 1 , ^ 

?l) " ~ ' ^r ( i- e) (^-Wi-ifci) 2 I (MP 1 • (427) 



This last integral cannot be expressed in terms of elementary functions at arbitrary e. In 
the limit e — > 0, instead, at fixed nonzero g it becomes 



,2\ / „I2 



/ + (?i^ 1 ) = -ln^Jlii^J . (4.28) 

Note that for I+(qi, q[) the limits e — > and q — > are noninterchangeable. However, it does 
not matter for J7~2, where I+(qi, q[) enters with the coefficient ~ q 2 . At e — > we obtain 

J2 = {ql + q? - q 2 ) (~ + g - C(2)) + y m(-g 2 ) + ?1 2 ln(-g-) - q 2 ln(-q 2 )) 

Appearance of the pole at e = in JT 2 does not contradict the statement that non-integrable 
infrared singularities in the limit e — > are absent in — x)]) + . The pole term comes 

from the ultraviolet region. Indeed, one can see from the expression (4.12) for A 2 that, after 
averaging over the azimuthal angle at D = 4 in the limit of large \k\\, we have 

I— ^ ~^( x (l-x)-2). (4.30) 



\x(l — x) ) + k\ 

The pole term in Eq. (4.29) is just the doubled result of the integration of the expression 
(4.30) over dxd 2+2e ki/n. Evidently, the ultraviolet divergency is artificial and appears to be 
the result of the separation of F s (ki,k 2 ) as shown in Eq. (4.9). Indeed, as it can be easily 
proved from formulas (4.12) and (4.13), the terms leading to such divergencies cancel in the 
sum A 2 + A 3 + (&4 <-> k 2 ). This means that the pole term in J 2 is cancelled by an analogous 
term in J% (see below). 
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4.3 Calculation of J% 



The integral J% is much more complicated than the preceding ones. As a consequence, the 
trick of integrating first over x, applied before, cannot be used in the calculation of 
because it leads to terms with the denominators containing a third power of k\. Such terms 
cannot be integrated over k\ by known methods. This complexity is connected with non- 
planarity of diagrams represented by J%, which is seen from the denominator i\t 2 related to 
the cross-box diagram. The complexity of contributions of the cross-box diagrams is well 
known since the calculation of the non- forward kernel for the QED Pomeron [17] which was 
found only in the form of two-dimensional integral. In QCD the situation is greatly worse 
because of the existence of cross-pentagon and cross-hexagon diagrams in addition to QED- 
type cross-box diagrams. It requires the use of additional Feynman parameters. At arbitrary 
D no integration over these parameters at all can be done in elementary functions. It occurs, 
however, that in the limit e — > the integration over additional Feynman parameters can be 
performed, so that the result can be written as two-dimensional integral, as well as in QED. 
Details of the calculation are given in Appendix D. The result is 

Js = ^(q 2 i + q? - q 2 ) Q + H-k 2 ) + l ) + J ^ w ?) + (ft <- ft') 



with 



J( qi , q 2 ; q) = [ dx f Q dz ((2 - XlX2 ) In (^j - 1 In ^ ) 

(x 2 qiq[ (ri(gi - r 2 )) - gf 2 gir 2 ) ^ 



+ (z(l - z)q' 2 2 q l q' l + q[ 2 (z qi k + (1 - z)q 1 q' 1 )) ^ 



Qo 



+^ (q?qi ( r i - 2gi) + 4xi^(^r 2 ) + q'tfiiq'tfi - ft: 7 "! ~ 91^2) + %i?"i)(9ir 2 ) - 2(q' 1 r 2 )(q 1 r 1 



+q{ 



12 



.(•■1 ■ j r\ \fJ%/ -''I 



' ' 2—(q 1 r 2 )q' 1 k + a; 2 (gir 2 )(g 2 - /c 2 ) + 2(g 2 r 2 )gig) + f — (qir )q[k 



1AQ 2 V xi 



+~T ( ~2 ln (%) ~ Tri I ( 2 ~ (ft^XVi + A;)r 2 - 2((x 2 gi + q 2 )r 2 )q x r 2 
?2 W V ^2/ Q J v x i 



1 ^i.^W 



>. " l 1 ■» I 771 ) (2 — (7: '(,)! 7' • /'•)'•() 



q\ 



d 



1 Q 2 
k 2 d 



1 



1 ^ 



- -f ( ( T5 + ^ ] - (q 2 r 2 )(q' 2 k) I - ^£ 1 - (q 2 k)(q' 2 n) I -j - ^£ 



(4.31) 



16 



Here 



n = zxq x + (1 - z)(xk - (1 - x)q' 2 ), r 2 = z((l - x)k - xq 2 ) + (1 - z)(l - x)q[; n+r 2 = k, 

g2 = _ x{1 _ x){g 2 z + g /2 (1 _ z)) _ z{l _ z){q 2 x + g /2 (1 _ x) _ q 2 x[l _ ^ £ = Q 2_ ^ 

r = zk + (1 - z)q[, Q 2 = -z(l - z)q' 2 2 , /j, 2 = -zk 2 - (1 - z)q' 2 , 
d = fi 2 fx 2 2 + k 2 Q 2 = z(l-z)x(l-x) [{k 2 - q\ - q' 2 )(k 2 - q' 2 - q\) + k 2 q 2 ) + qfq 2 xz(x + z - 1) 

+^ 2 (l-x)(l-,)(l-x-,), £ = ln(-^L). (4.32) 

5 Non- forward kernel 

In pure gluodynamics, which is considered here, the part K.^ of the BFKL kernel (2.2), 
related to the production of real particles, for any representation 1Z is given by the sum of one- 
gluon K G and two-gluon IC GG contributions. Using for the last of them the decomposition 
(3.46) we have 

Kf > = /Cg*> + A(a R - b R )K% G + bnKol , (5.1) 

where the colour group coefficients a R and b R are defined in Eqs. (3.25) and the one-gluon 
contribution Kffi is given by Eq. (2.12) (for arbitrary D see Eq. (4.10) in the last of 
Refs. [8]). The two-gluon contribution for the octet channel IC GG was calculated in the 
second of Refs. [10] (see there Eqs. (61) and (63) for arbitrary D and for D — 4, respec- 
tively). The calculation of the "symmetric" contribution IC GG performed in this paper solves 
the problem of finding the expression of the non-forward BFKL kernel for all possible colour 
states in the t-channel. This contribution is determined by Eq. (4.14), where Jq — and J i: 
for i = 1 + 3, are given by Eqs. (4.21), (4.26), (D.8) and (4.22), (4.29), (4.31) for arbitrary 
D and for D = 4 correspondingly. Note that everywhere in these formulas the bare coupling 
constant g is used. The transition to the renormalized coupling in the MS scheme takes 
place by means of Eq. (2.5). 

For the most important colour singlet case, using C\ = a x = 1 and c 8 = b± — 1/2, from 
Eqs. (5.1) and (2.12) we obtain 

^ = + . (5.2) 

Because of the significance of this case let us consider it in more detail. We present the 
kernel K,^ obtained from the above stated sources in the limit D = 4 + 2e — > 4 as sum of 
two parts, 

£(1) = x-sing + fcreg ^ ^3) 

where the first, given by 

K g (g^ = ^r(i- e) ( p q )\ 1 + 9 » 
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3e 



k2 \ f 11 67 n ,/M / 404 11 

i?J {-^ + ¥- 2 « 2 ) + £ (-i7 + 14 « 3 » + T 

^2 _ „, /\,2> 



(5.4) 



contains all singularities and the second, putting e = and g 2 = a s (fi )N c /(4ir) , becomes 

2 (9 - + g-g)(c(2)-f)-H( ??lD (g 




-(» ta (R + ^(8)-(S)(? ta (fi-f-(s: 

q?% 2 + qiQl 2 ^ fg\ + ff 2 (ff - 3g«) ln (^] + ^ W " # 




2A; 2 V^ 2 / 2fc 2 \? 2 2 y 2A; 2 

i i ^2/T2 ^2 ->2\ , o^2^2 (q 2 ~~ q 2 )(q 2 q2)(qi 2 ~ q^ 2 ) , ^2^/2 , ^2^/2 
+ 5 {k - q x - <? 2 ) + 2gi g 2 ttt + ?i ?i + 9 2 9 2 



- y (# 2 + «£ 2 ) ] ft 2 , ft 2 ) - 2 J(ft, &; 9) - 2J(-£, -ft) 



2fc : 

'+ (ft <—>#!. (5-5) 



Here the functions J(/c 2 ,g|,g 2 ) and J(qi,q2',q) are defined in Eqs. (4.23) and (4.31) corre- 
spondingly. 

All singularities of A^ 1 ** are contained in the first part. We remind that Kffi and k!q G 
separately contain first and second order poles at e = (see Eq. (2.12)). In the sum of these 
contributions the pole terms cancel, so that at fixed nonzero k 2 , when the term (k 2 / jj 2 ^ in 
Eq. (5.4) can be expanded in e, the sum is finite at e = 0. But the kernel (5.4) is singular 
at k 2 = so that, when it is integrated over q 2 , the region of so small k 2 values such that 
e| In (k 2 / n 2 ^ | ~ 1 does contribute. Therefore the expansion of (k 2 / '/j 2 ^ is not done in 
Eq. (5.4). Moreover, the terms ~ e are taken into account in the coefficient of the expression 
divergent at k 2 = in order to save all contributions non- vanishing in the limit e — > after 
the integration. 

As it was already discussed, the "symmetric" part K.%% of the kernel (5.1) is finite in 
the limit e = 0. Moreover, it does not give singularities at e = when the kernel is used 
in the equation for the Green's function. Indeed, the points q 2 = and q 2 ' = do not give 
such singularities due to the "gauge invariance" properties (2.13), because these properties 

(s) 

are valid for any representation. Its validity for fC GG can be checked explicitly using the 
properties of b a/3 . It is not difficult also to see that fc GG has not non-integrable singularities 

— * 

in the limit e = at k = 0. 

For the singlet case the infrared singularities of /C^ 1 ** must be cancelled by the singularities 
of the gluon trajectory after integration of the total kernel with any nonsingular at k = 
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function. The total BFKL kernel in the singlet case must be free from singularities. It is not 
difficult to see that it is the case, using the equality 



^H-^(i + mg))- 5 ; 



11/1 . 2 f-t\\ {67 



(5.6) 



It is convenient to represent the total kernel in such a form that the cancellation of singu- 
larities between real and virtual contributions becomes evident. For this purpose let us first 
switch from the dimensional regularization to the cut-off k 2 > A 2 , A — > 0, which is more 
convenient for practical purposes. With such regularization we can pass to the limit e — > 
in the real part of the kernel, so that its singular part assumes the form 



lC S r in9 (qi,q2-,q) - - <f 2 ) 2 - A 2 )/C^(fi, <? 2 ; q)\ £=0 



a s (fi 2 )N c f q?qp + q( 2 q 2 2 



2tt 2 



k' 



x 1- 



a s (n)N c III Ik 



67 



^lT ln (^]-T +2C(2) )r (( *-« A) - 



(57) 



The trajectory must be transformed in such a way that the cut-off regularization gives the 
same result as the e regularization: 



a 2 \ 



In 



A 2 



« s (/i 2 )iV c 



47T 



11 



In 2 



-t 



-In 2 



'67 7T 



3 ; ln vp 



6 \ V A 4 

+ 6C(3) 



(5.8) 



v/i 2 ; / v 9 

It is easy to check that the integral over d 2 q 2 of the total kernel (2.2) with any function 

— * 

non-singular at k = is A -independent in the limit A — > 0. Moreover, it is equally easy to 
find a form of the kernel which does not contain A . It is sufficient to find a representation 



u x {-q 2 ) = J d 2 q 2 Uqi,q2)0{{qi ~ ^ ~ A 2 ) 



(5.9) 



with such a function f w that the singularity non-integrable at k — qi — q 2 — q{ — q 2 = is 
cancelled in the "regularized virtual kernel" 



K r v 9 {qii q.2\q) = 92) + /«(&', qi) + 

After that we can proceed to the limit A = 0, obtaining 



JC s r ing (qi,q2-,q)\e=o 



— *9 — */ 9 
9 2 9 2 



(5.10) 
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Of course, the choice of the function f u contains a large arbitrariness. One simple choice is 
a s (fi)N c ( 11 /A; 2 \ 67 / H A / A fc! 



x 1- 



4, vH^H +2C(2)+ ( 6C(3) ~>)^ r (512) 



We have to say that the integral (4.31) for J(qi,q2',Q) entering into all kernels besides 
the octet ones definitely is presented not in the best form. We have decided to present it 
in such shape in order to give a possibility of further development to people interested in 
this subject. The results of our efforts on simplification of the kernel and investigation of its 
properties will be presented in a subsequent paper. 
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Nazionale di Fisica Nucleare - gruppo collegato di Cosenza for their warm hospitality while 
a part of this work was done. 



Appendix A 



For the colour group SU (N c ) with N c = 3 the possible representations 1Z are 1, 8^, 8s, 10, 10, 27. 
Corresponding projection operators are 



(bb'\Vi\ad 



a 8 bb >5 aa / 



Nl - 1 ' 



{bb \V 8a \aa ) = 



(bb'\V 8s \aa') = d Wc d a 



N r . 



''N? -4 ' 



(B'lAoK) 

(bb'\VjQ\aa') 



* to <W - 6 ba ,5 b , a - -f Wcfaa , c + if ba , c d Vac + id„ e f Vl 

hah'a 1 — ha'h'a — J-j~ fbb'cfaa'c — ifba'cdb'ac — id ba 'cfb', 



{bb'\V 27 \aa'} = ^ 



2 \ 2{N C + 2) 

1 + ~rr ) (Sbadb'a' + $ba'$b'a) —Shh'Sn 



N C (N C +1) 



(Al) 

(A.2) 
(A.3) 
(A.4) 
(A.5) 
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1 + 



N c + 2 



dbb'cdaa'c + dbacdb'a'c + db'acdba', 



Here f a bc and d a bc are defined by the relation 



t a t b = (d abc + if abc y- + 5, 



ab 



2Nr 



(A.6) 



(A.7) 



where t a are the group generators in the fundamental representation, normalized by the 
requirement tr{t a t h ) = S a b/2 and / is the identity matrix. 

For generality, we do not put here N c = 3, so that above expressions are valid for the 
SU (N c ) group with arbitrary N c . Corresponding representations in this case have dimensions 



rai = 1 , n 8a = n 8s = Nj-l , n 10 = 



(N?-4)(N*-1) 



n 2 7 



However, at N c > 3 there is an additional representation with dimension 

(iV c + l)iV c 2 (jV c -3) 



™iV c >3 



(N C + 3)N^(N C -1) 
4 

(A.8) 
(A.9) 



and projection operator 



(bb'\V Nc>3 \aa') = - 



+ 1- 



1 2 \tx A 4- A A n , 2(N e -2) , A 

J- — ~rr ) [ObaOb'a' + Oba'Ob'a) + T7TT7 —°bb>O a a> 



N C (N C - 1) 



dbb'cdaa'c dbacdb'a'c dba'cdb'ac 



In SU (3) this projection operator turns into zero due to the equality 

dbb'cdaa'c + dba'cdb'ac + d bac d a 'b'c = g (<W<W + ha'h'a + ^6o^a'6') 



which holds at iV c = 3. The following useful identities with 



T bc ~~ */a6o -Dfcc — rfabc, 

are valid at arbitrary N c : 



pa pb 



"ifabcF t 



pa D b 



ifabcD C 



(A.10) 



(A.ll) 



(A.12) 



Tr (T a ) = Tr (D a ) = Tr (T a D h ) = 0, Tr (r a T b ) = N c 5 a \ Tr (D a D h ) = ^ — -5 ab , 



N, 



Tr (T«T b T c ) = t'-yfabc, Tr (T a T b D c ) = ^d abc , 



Tr (D a D b T c ) 



N 2 C -4 
' 2N r 



fabc, Tr (D a D b D c ) 



N? - 12 



2N r 



'd a bci 



Tr (r a T b T c T d ) = 5 ad 5 b c + \ (fofa + S ac S bd ) + ^ (f adt f ba + d adl d b ci) , 



Tr (r«T b T c D d ) = & (d adi f bci - fasAc^ , 



21 



Tr (T a T b D c D d ) = X - (5 ab 5 cd - 5 ac 5 bd ) + ^l^f adi f bci + ^d adi d bci , 
Tr (T a D b T c D d ) = - X - (5 ab 5 cd - 5 ac 5 M ) + ^ (jadifba + d adi d bci ) , 
Tr (T a D b D c D d ) = i^f adi d bci + i^-^f abi d cdl + i^d abl f cdl , 
Tr (D a D b D c D d ) = -^L_8 ad 5 bc + -5 ac 5 bd + -^__5 ab 5 cd 

~r ^ Jadijbci ~r ^.jy ^adi^bci ^abi^cdij 

2 

jadijbd + d adi d b d — jabijcdi — dahidcdi + — (5 ad 5 bc — 5 ab 5 cd ) = 0. (A. 13) 

These equalities (A. 13) can be derived from the relation (A. 7) and the completeness of the 
matrices t a and /. The properties of the projection operators 

ViVj = 5^, E^ = / (A.14) 

i 

can be easily obtained with the help of these equalities, as well as the coefficients c-r: 

1 1 /AN 

ci = 1 , c 8a = c 8a = - , cio = cjo = , c 27 = -c Nc>3 = -— . (A.15) 



Appendix B 

Here and below all vectors are taken transverse (D — 2) -dimensional, although the transver- 
sality sign _L is omitted. If the vector sign is not used, the Minkowski metric is assumed, so 
that (ab) = —ab. We use a standard representation of logarithms, i.e. 

hia=4-a u \ v=0 (B.l) 
dv 

and the Feynman parametrization 



n o = » 1 y n / n ***r 1 — ^v- ■ (b-2) 

nrK) w ° y /„ N E< 



a,- 



i=l ^E J i=1 



Using the notations 

R = {k 2 - yk x ) 2 -y(a x - yb x ) , k x = k(l - x) - q 2 x = q x (l - x) - q 2 = k - q x x , 



/2 

= - (k 2 (l -x) + qjxj , 6 S = -k 2 £ = a x + qfx(l - x) , Ji = -y Ji + <-> g-) , (B.3) 
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from Eq. (4.20) we obtain 

d f d 2+2e h ((qi-h) 2 ^ 

d r 1 



q\q\ 



dvJ 7^+^(1 -e) (klYk 2 (qi - h) 2 



g| + 2kiq 2 \ 
k\ ) 



v=0 



du Jo ^ Jo 



dx{l-x) v x- u 



-I- 

I/) J 7T 



d 2 ^h (2yq 2 q 2 



q\v 



r(l + v)Y{\ - u) J 7r 1+e r(l - e ) v R 3 (l-x)R 2 

d_ r 1 dx{l-xfx- u r 1 dyy'- 1 
~Yvk r(l + v)Y{\ - v) Jo (ax-ybx) 1 -' 



+ 



(qj - k 2 - 2k 2 q 2 )v s 
xR 2 



x 



(1 - e)q\q\ vq 2 



x 



(1 ~ e)qlql 
(a x - yb x ) ~ £ 



In 



a x - yb x ) (1 - x) x 

= f dx [ dyy 6 ' 1 
Jo Jo 

d q\ 



+ -{q\ - k 2 - 2yk x q 2 ) 



\v=0 



X 

l-x 



ln(l — x) 



dx (a x - yb x ) 



lnx 



d q\-k 2 - 2y k x q 2 
dx (a x - yb x ) 



l-e 



(B.4) 

Subsequent integrations can be performed only in the limit e — > 0. Note that proceeding 
to the limit e — > in the integrand leads to a wrong, although convergent, integral. The 
singularity at y = requires an accurate consideration. After integration over y we obtain 



Ji 



dx 



q\q\ 



In 



1 — x 



In 



^x &x 



+ 



In(l-x)^ln 

dx a r 



&x b x 



lnx 



d_ ( q\ - k 2 _ 2k x q 2 
dx 1 a T b T 



In 



Q>x &x 



X 



1-x 



+ q 2 (k 2 - q 2 2 ) ln(l - x) + (tf - r)(A; 2 - g 2 ) lux In a* 



(B.5) 



Last terms come from the singularity at y — 0. Performing appropriate integrations by 
parts, after simple though slightly tedious steps we arrive at 



Jl = Jo dx (i= 



'^ 2 + lW- 



/c 2 



+ 



1 



Using the equalities 

xdx 



(<£(<£ - ql - e)x + e(gl + ql) - (ql + In (^-)) ■ (B.6) 



Jo 



x ln 



(fc(l — x) — g2^) 2 

k 2 (l — x) + g 2 x^ 
g 2 x(l — x) y 



ln 



fc 2 (1 — x) + g 2 2 x\ q 2 + k 2 — q 2 f 1 



q 2 x(l — x) 



2?i 2 



/ 
Jo 



dx 



(k(l — x) — g2^) s 



1 / t 2 « 2 

^(l-^)-^(l-^) 



2gi 



92 



i in O in (^ 



L(x) = -Li 2 (x) = / -^ln(l-y), Li 2 (l - x) + Li 2 (l --) = -- In 2 x , 
Joy x l 



(B.7) 
(B.8) 



i/=0 
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we obtain 



2A; 2 Jo K, \a T 



+ ^J^(5WSb (B,) 



M. 



This result can be reached by another way, using the analyticity properties of J\. To do 
that, let us present J\ at e = as an integral in the Minkowski space: 



Ji 



dz 



d 2 h 



1 



1 



wr(fc| + iO) \z-k 2 -i0 z- ( qi - k^ 2 - iO j 

q\ + 1k x q2 



x ^2 ^2 H2 T ^\H2 (Rim 

^(^ + io)(( 9l -fc 1 ) 2 + <o) ((q.-hy + to) (kf + io))- 1 • UJ 

Here fci, q\ and g 2 are considered as vectors in the two-dimensional Minkowski space, so that 
we have 

d 2 h = dk^dkP, k 2 = (k^) 2 - (k^) 2 , (B.ll) 

and so on. Eq. (B.10) determines J\ as a function of q 2 , q\ and k 2 for arbitrary values of 
these variables. For q\ = —q 2 < 0, q 2 = —q 2 < and k 2 = —k 2 < Eq. (B.10) turns into 
the function given by Eq. (4.20), that can be easily seen by making the Wick rotation of 
the contour of integration over kf" 1 and performing integration over z. At fixed negative q\ 2 
Eq. (B.10) determines the real analytical function of k 2 with the cut at k 2 > 0. According 
to the Cutkosky rules, one can find a discontinuity on the cut rewriting Eq. (B.10) as 



Ji 



d 2 h 



m{kl + iti) \\{{qi -h) 2 + i0) 



q\ - 2k x q 2 



1 



1 



+ 



z \k% + i0 k 2 -z + i0) 
1 1 \\ 



kl + iO (q 1 - h) 2 -z + iOJ ((gi - h) 2 +i0) \kf - z + iO ( qi - h) 2 -z + i0 / 



omitting the last term and making the substitutions (assuming k^ > 0) 



(B.12) 



ki + iO kl + iO 

1 1 

k 2 2 + i0k? - z + iO 



(-2m) 2 5(k 2 )5(k 2 )6(k^)6(kr) , 



'-2m) 2 5{k 2 )5{z - k 2 )9(k ( 2 0) )9(k ( ">) 



(B.13) 



Using these rules and removing the 5-functions by integration over k\ (the most appropriate 
system for this is k^ = 0, k 2 = (k^) 2 ), we obtain for the imaginary part 



roo 

s j\ = it j dz 

1 k 2 6(k 2 



i=± 



J_ k 2 9(k 2 -z) \ k 2 -q\-q\ 

' n /to \ T 



k 2 yn^n® — z) k®z z(k 2 — z)ni 



+ 



+ 



+ 



+ 



6{k 2 



k 2 



K?-Z 



z(k 2 — z) ) k 2 ^ — z) k 2 z k 2 (K® — z) [k 2 — z) \k 



Z 



(B.14) 



24 



where 4 and k± are given by values of (qi — ki) 2 on the mass shells k\ — , k\ = and 
= , k\ — z respectively, so that 



K± = 



+ ql -k 2 ± ^(ql + ql - k 2 ) 2 - 4qlql\ , «± = 4 + ^ - 4) • (B.15) 



The integration over z is quite elementary and gives 

'k 2 - q\ - q\ + yltW-ql-qlY-Witi 



7T 
2k 2 



(fc 2 - gf - ql) 2 - Aq\ql In 



k 2 - ql - q\ - ^k 2 -q 2 -q 2 2 ) 2 -Aq 2 q 2 , 



+ (k 2 - q\ + ql) In 



5? 



(B.16) 



The use of this equation and the equality (see Refs. [10] and [15] ) 

'k 2 - ql - q 2 + ^V-ql-qlY-Aqlql' 



71 J (k 2 - ql- ql) 2 - Aqjql In 



3 Uqfql - (k 2 - ql - q 2 2 ) 2 } f 



k 2 -ql-ql- J (k 2 - ql - ql) 2 - 4qfq% t 



dx 



In 



' k 2 (l — x) + q 2 x~ 



(B.17) 



lo (k(l — x) — q 2 x) 2 \ q 2 x(l — x) 

gives the result (B.9). The absence of a polynomial in k 2 can be easily checked by considering 
the integral (4.20) at large k 2 3> q 2 . Integration regions which could contribute in this case 
are ~ \k — k\\ ~ \k\ and \k\\ ~ \q\ — k\\ ~ \q\\. But the first region gives a vanishing 
contribution because of the smallness of ln((<fi — ki) 2 /k 2 ) there. In the second region the 
integrand of Eq. (B.9) with the required accuracy is anti-symmetric with respect to the 
exchange k\ <-> qi — ki, so that its contribution vanishes as well. 



Appendix C 



We start from the integrals 

\ L+ c + d L. 



dx 



o \x(l — x) J led led b — V Jo \ x(l — x) 



' 4 W 



2W 



L b + b ' L 
L+ ~b-b' L ~ 



1 j / x 2 (x 1 tit' 1 
dx ' 



i ^ — i 1 



dx 



x(l — x) 
4(iit[)- v 



a(b - b') 



' So 



1 dx ^ ^^(^i^i^i) 



x(l — x) 



1 

2^2 



L C + C L 
~ L+ ~b-b' L ' 



lo \ x(l — x) 
where 



1 

W 



. 1 + |± (W+a(6+60) __^_ 



, a(b-b') 2 . 



a — k 1 , c — (q± — ki) , c' — (q[ — k±) , b = c — a, b' = d — a , 



(C.l) 



(C.2) 
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and after some algebra we arrive at Eq. (4.25). We remind that in the integrands we omit 
the terms giving zero after the subsequent integration. Then we use the following equalities: 

f d 2+2e ki a fc\ / 2 y+i T(l + e)r(2 + e) 
J n i+*r(l - e) b \a) ~ V 91 J er(4 + 2e) ' 

f g^fci i , /c\_ r 2y r 2 (i + 6) 

7 7ri+T(l-e)& W ^ ^ er(2 + 2e) ' 
/■ rf 2+2e A;i (gifci) , fc\ t 2 u r 2 (l + e), ((i , ^ xx 

/ ^w^) V ln U) ^ H) ^tt4 m + e) - ^ (l + 2e)) ■ 

,2/ ^ 2+2£ fei (gi-Wi-fci) ln /c\_ , / 2 y r 2 (i + e ) 

?l i ^(l - e) ac ln U ~ HO e 2 T(l + 2e) ' 

2 f d 2+2e fci L_ / 2 y+i T 2 (l + e) 

9 7 ^f(l - e) b - V ~ H J e r(2 + 2e) ' 
r d 2+2e fci c + c' / 2 y+i T 2 (l + e) 

J 7r^T(l -e)b-V ~ = V g J e r(2 + 2e)r(3 + 2e) ' 

The first three of these integrals can be easily calculated with the help of the representation 

'tofi^U/ 1 -^-, (C.4) 
o \ a y jo a + ox 

the forth with the help of Eq. (B.l) and the last three using the representation 

^-=r d 4 -. (c.5) 

b-b' Jo cx + d(l-x) v ; 

Using these integrals we arrive at Eq. (4.26). As for the integral I+(qi,q[) of Eq. (4.27), it 
can be written as 

InfM). (C . 6) 



7r 1+e r(l - e) (51 - /ci) 2 (gi - ki) 2 \q 
The first integral in Eq. (C.6) can be easily calculated at arbitrarye; we find 

d 2+2e h f( qi - hM - h) ^ / cd \ c + c'^ fkT 



1+e r(l-e)V cd \(q 2 ) 2 J cd \q 2 ) 

= j^+ae) ( 2 ("^ (l " + ^ (1 " e) " ^ (1 + e) + ^ (1 + 2e l 

4(K 2 ) e+ H 2 ) e ))- ^ 

The second integral in Eq. (C.6) was analyzed in Ref. [16]. In the limit e — > we have 

/■ rf 2 ^ g 2 /A; 2 \ _ 1 / 2 y /(g 2 ) 2 \ 1 2 / g 2 \ 

i vr^r(l- e) M ~ 7 H J 111 fej - 2 M IgfJ • (C - 8) 
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Appendix D 



We use Eqs. (4.15), (4.13) and after a simple algebra we obtain 



7 = f 1 dx( 1 f d2+2eh 1 
^ 3 7o V^(l - x) J 7^+^(1 - e) t~i> 



(D-2) 



x ia ; 2 (^- 2^0(^-2^ fc 2 ) 



_:C2 -(?i?i)(^i(?i - fo)) + gf (2^-^ - 2^ + gi/ci ) + Axiql(q[k 2 ) + q[qi(q[qi - q[h - q x k 2 ) 



X\ \ X\ 



+2(q[k 1 )(q 1 k 2 ) - 2(q' 1 k 2 )(q 1 k 1 ) + || (2^(91*2) (q[ h - kk 2 ) 



+x 2 (q[k 2 )(ql - k 2 + 2q 1 k 2 ) + 2(q 2 k 2 )(q l q + gifc 2 )) + gf g 2 



2 (^k ^ig^h] 
21,2 



(i).n 

Terms in Eq. (D.l) with x\ in the denominators require subtraction, so that the prescription 
(3.41) is important for them (and only for them). Note that neither the integrand in Eq. (D.l) 
itself, nor the subtraction terms contain non-integrable infrared singularities. Nevertheless, 
the integral J% has a pole at e = 0. The pole comes from the ultraviolet region. Evidently, 
the ultraviolet divergency is artificial and appears as a result of the separation shown in 
Eq. (4.9). It is easy to see from Eqs. (4.12) and (4.13) that the terms leading to such 
divergencies cancel in the sum A 2 + A 3 . 

The integration over k\ is performed using a standard Feynman parameterization. The 
basic integrals are 

r d 2+2e h 1 r 1 dz r d 2+2€ h k[ r l dzr? 

J ^Y{l-e)W 2 ~ XlX2 Jo Q 2 ^ ' J 7r 1+e r(l -e)tj 2 ~ XlX2 Jo Q 2 ^ ' 

~ x i x 2 ^ 2 n_ f) [nr 2 - — Q , (D.2) 



7r i+e r(1 _ e) ij 2 — y Q2 ( i- e) v 1 2 26 

where 

T\ = zxqi + {l—z){xk—{l—x)q 2 ) , r 2 = z((l—x)k—xq 2 )+(l—z)(l—x)q 1 ] r\+r 2 = k . (D.3) 
The integrals with t\t' 2 k 2 can be calculated joining first t\ and i 2 and then kf. We have 

/• d 2+2e h 1 _ _ r 1 f 1 ydy 

J ^Y{l-e)ht' 2 k 2 ~ XlX2[l €) J dZ Jo (y(tf + yr 2 )) 2 - ' 

f d 2 ^h k[ f 1 d yy 2 < 

J ^(1 - e) U' 2 k 2 ~ XlX2[L €) Jo dZ Jo (y^ + yr 2 )) 2 -* 

r 1 dzr? 

'0 1W 
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d 2+2e h fcffi 

= —x\%2 \ dz 
Jo 



£12:2(1 -e) f dz f — — 
Jo Jo (via 



ydy 



(y (jj>i+yn)) 



^2^7 \ v n r i - 



v g^y + yr- 



f rr / 1 



In 



7 /ii/ 



'91 



In 



2(1-6) ; 

+ 0(e) , (D.4) 



where 



A*i = <5 2 - , lA = -zxq\ - (1 - 2;)(x/c 2 + (1 - x)g 2 2 ) , 

/1 2 = -z((l - x)k 2 + xq 2 2 ) -(l-z)(l- x)q[ 2 . (D.5) 
Finally, the integral with t\t' 2 k\k\ can be calculated joining t\ and t 2 , then the result with 



k\ and subsequently with k\. We obtain 

d 2+2f h k^k\ , n W1 , f 1 f 1 f 

^w^wm =XlX2{ e){ ] Jo 1 Jo 



dzydyt 3 dt 



(t(y(ri + 2kr 1 )-k 2 + t(k-yr 1 ) 2 )) 



3-e 



X 



(k - ynYik - t(k - yr^Y + -f— - (y (/i 2 + 2fcn) - A; 2 + - yrO 2 ) 



1 dz 
= xix 2 I — 
10 d 



Q 



t4 



k " kU 77 + ^ " ^k v \- 2 ~ ^-C - fc^r? ^ - ^£ 



A; 2 d 



d 



k 2 



^ [Q 2+ ^ + 2 C 



+ 0(e), 



(D.6) 



where 



d = fi 2 lf i 2 2 + k 2 Q 2 = z{l-z)x(l-x) ((A; 2 - q\ - g 2 2 )(A; 2 - gf - q\) + A; 2 g 2 ) + g 2 g 2 a;z(x + 2-l) 

(D.7) 



+q[ 2 q 2 2 (l -x)(l -z)(l -x-z), C = In 
At arbitrary D we have 

J3 — [ dx [ dz 
Jo Jo 



-k 2 Q 2 



q iq [ [x lX2 Q 2 ^ + -(x 2 Q 2e - Qo)) 

(x 2 q 1 q' 1 (n(gi - r 2 )) - gf gi r 2 ) 
+ (z(l - z)q 2 q x q\ + gf {zq x k + (1 - z)gigi)) 2 (i- e ) + Q2(i- e ) (^i ( ri ~ 



^Q^ej x i x M - 2giri)(gf - 2q[ r 2 ) - — 



+4rig 2 (^r 2 ) + q[qi(q[qi - q[n - gir 2 ) + 2(q[r 1 )(q l r 2 ) - 2(q' l r 2 )(q l r l )^j 



+Q? I ydy - - !_ e (-— qi{q[ + k) + x 2 qiq[ + gig 2 ) - 



1/(1 -e) 



.^2 



x 2— (gir 2 )(^fc - y(gj + fc)r 2 ) + £2(<7i?" 2 )(g 2 - k + 2yq l r 2 ) + 2(g 2 r 2 )(gi(g + yr 2 )) 

Xi 



+ 



(y^l + yri)) 6 Vzi 



1/(1 " e) 



(y(^ + yr 2 )) 2 e Vx! 
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(qiro)(q'ik - y(q[ + k)r ) 



+(2-e)(l-e)g? / 
Jo 



o (t(y(i4 + 2Atx) - k 2 + t(k - yn) 2 )) 



[(q 2 (k - yr l )){q' 2 {k - t(k - yn))) 




yn) 2 ) 



I 



+ (ft <-»• ?D , 



(D.8) 



where 



Qjj = - fx 2 = -zk 2 -(l-z)q[ 2 , ro = zk + (l-z)q' 1 , r 2 = Q 2 - fi 2 . (D.9) 

In the limit e — > we arrive at Eq. (4.31). 
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